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Ultra-clean crystals of Sr3Ru207 undergo a metamagnetic transition at low temperatures. This 
transition shows a strong anisotropy in the applied field direction with the critical field He rang- 
ing from ~ 5. IT for H perpendicular to c to ~ 8T for H || c. In addition, studies on ultra-pure 
samples revealed a bifurcation of the metamagnetic line for fields in c-direction and it has been 
argued that a nematic phase emerges between the magnetization jumps. The aim of this study is 
to explain the field-direction anisotropy of these phenomena. Based on a microscopic tight-binding 
model, we introduce the metamagnetic transition by means of a van Hove singularity scenario. 
We show that the rotation of the O-octahedra around the c-axis observed in this material intro- 
duces a staggered spin-orbit coupling within the planes and naturally leads to an anisotropy in the 
low-temperature behavior around the metamagnetic transition. In particular, the low-temperature 
(nematic) phase is affected. We show that uniform in-plane magnetic fields induce a (commensu- 
rate) staggered magnetic moment component which can suppress the low-temperature phase. In 
contrast, the response to fields along the c-axis remains unaffected and thus, also the corresponding 
low-temperature phase . As a concrete example, we choose a nematic Pomeranchuk instability for 
the low-temperature phase. An experimentally testable prediction of this work is the occurrence of 
a staggered magnetic moment in response to a uniform magnetic field perpendicular to the c-axis, 
which should be accessible by neutron scattering. 

PACS numbers: 71.10.Hf, 73.22.Gk, 74.70.Pq 
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I. INTRODUCTION 

The ruthenium compounds of the Ruddlesdcn-Popper 
series, Sr„+iRu„03,i+i, have been the subject of inten- 
sive research for over a decade due to their interesting 
ground-states. While the quasi-two-dimensional single- 
layer (n=l) is an unconventional superconductor likely 
with p-wave pairingj^^^ there is a growing tendency to- 
wards ferromagnetism with increasing layer number. The 
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FIG. 1: Schematic phase diagram of Sr3Ru207 for fields H 
applied with an angle 6 versus the afc-plane. The surface 
represents first-order transitions separating a region with low 
(spin) polarization from a region with high polarization. The 
thick black line connecting {9,H,T) = (0°, 5T, 1.2A') and 
(80°,8r, OA") is a line of critical endpoints. For details see 
main text. 



three-dimensional infinite-layer compound, SrRuOa, in- 
deed realizes an itinerant ferromagnet.i^The bilayer com- 
pound (n=2) with its intermediate dimensionality, how- 
ever, shows no ordering down to lowest temperature^ 
and is supposed to be a strongly correlated Fermi liquid. 
Still, applying uniaxial pressure can induce a ferromag- 
netic transition and this compound is, thus, expected to 
be on the verge to a ferromagnetic instability.SI This is 
also supported by inelastic neutron scattering- and band 
structure calculations.^ 

As was first discussed by Wohlfarth and Rhodes,'^ such 
proximity to ferromagnetism can result in metamagnetic 
behavior, a superlinear rise in the magnetization over a 
narrow region of applied magnetic field H. This phe- 
nomenon was observed in a number of system^ and also 
in Sr3Ru207 it was founcF'with a critical field He, how- 
ever, that depends strongly on the angle of the field ver- 
sus the a&-plane. While for 6* = 0° the critical field 
He ^ 5. IT, it rises to He ^ 8T for 6 = 90° . In ad- 
dition, a first-order transition occurs for in-plane fields 
below T* sa 1.25K, while there is only a crossover for 
fields parallel to c. It was therefore suggested that the 
field angle could be used as a tuning parameter for a 
line of first-order transitions that goes to T* = around 
9 = 80°, thus, realizing a quantum critical endpoinlP^ 
(see schematic phase diagram in Fig. [T]). 
However, when trying to reach this quantum critical end- 
point on ultra-pure single crystals with residual resistiv- 
ities down to pq < l/ii7cm, a splitting of the metam- 
agnetic transition into two jumps was observed. These 
jumps define an intermediate phase whose exact bound- 
aries could be determined by measuring of several ther- 
modynamic properties.'ii' Later, it was shown that this 
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phase breaks the symmetry of the crystaP' and it was 
argued that this was due to an induced anisotropic elec- 
tronic state with a symmetry-breaking Fermi surface de- 
formation similar to a Pomeranchuck instabilityP^ That 
this kind of a phase, also called nematic phase in anal- 
ogy to liquid crystal phases, can lead to two consecutive 
metamagnetic transitions had already been shown in a 
paper by Kee and Kim.l^ In addition, the anomalous T 
dependence of the susceptibility x and the specific heat 
coefficient 7 could be explainedl^Sl Moreover, the two-fold 
degeneracy of the nematic phase allows for domain for- 
mation, such that domain-wall scattering could account 
for the increased resistivity of the intermediate phasel^^ 
The electronic structure of a single layer of Sr3Ru2 07 
is dominated by bands originating from the 4d t2g or- 
bitals dyz^ dzx and dxy hybridizing with the O 2p or- 
bitals. This leads in a simple approximation to two 
quasi-one-dimensional bands with mainly dyz and d^x 
character and a two-dimensional band stemming from 
the dxy orbital. These three bands are then addition- 
ally split due to the interlayer coupling resulting in 6 
bands. An important consequence of the bilayer split- 
ting is that one of the two bands coming from the Axy 
orbitals is shifted closer to the van Hove singularity. This 
was also confirmed by recent angle-resolved photoemis- 
sion spectroscopy (ARPES) measurements.'^ A chemical 
potential in the vicinity of a van Hove singularity is the 
condition for the scenario described by Binz and Sigrisll^ 
for a metamagnetic transition. Proximity to a van Hove 
singularity can also lead to a nematic phase accompany- 
inga metamagnetic transition as described by Grigera et 
g^jjiT] The anisotropy in the critical field strength could 
then be explained by spin-orbit coupling (SOC) effects 
similar tcP^ leading to an anisotropic effective g-factor. 
A different route is taken by Raghu et al.!^ and Lee et 
who studied a model where the metamagnetic tran- 
sition comes from the two (bilayer-split) one-dimensional 
bands. The anisotropy is then again introduced by con- 
sidering spin-orbit interaction on the Ru-sites and the 
nematic phase can be understood as an orbital ordering 
among these one-dimensional bands. 
Both these routes suffer, however, from a short-coming: 
even though they can describe the existence of a nematic 
phase and a dependence of He on the field angle, they 
cannot explain why the nematic phase occurs only for 
fields almost parallel to the crystalline c-axis. 
We will address this point in the present work. We study 
a model based on a two-dimensional band in a single layer 
originating from the dxy orbitals. The bilayer effects are 
only taken into account by placing this band closer to the 
van Hove singularity. Starting from this, we will consider 
the effect of a lattice distortion in the planes. The O- 
octahedra in Sr3Ru207 are rotated by e.S'^^Hand we will 
show how this introduces a staggered spin-orbit coupling, 
an effect similar to the Dzyaloshinski-Morya interaction 
for localized spins, here, however, for itinerant electrons. 
For magnetic fields applied in the plane, this will add 
a component with wave vector Q = (tt, tt) to the static 
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FIG. 2: Different hoppings in the three-band model of Ru 4 
Axy and O 2 pa, / py orbitals. The relative signs of the hopping 
integrals come from the phases of the orbital wave-functions. 



susceptibility. The induced commensurate spin-density 
wave (SDW) will open gaps in the Fermi surface close 
to the van Hove points which will have an impact on 
the occurrence of any instability that emerges due to the 
proximity to a van Hove singularity. We choose here the 
electronic nematic phase to examine this aspect, since it 
relies on the presence of the van Hove singularity and 
represents one of the most promising candidates for the 
intermediate phase. In order to discuss the essential influ- 
ence of the spin-orbit coupling on the phase diagram we 
adopt here a mean-field approach, with the short-coming 
that critical fluctuations are not included well. While, 
in particular, quantum critical fluctuations represent an 
intriguing part of the phenomenology of this metamag- 
netic transition, we assume that they are not essential to 
understand the basic effects due to spin-orbit coupling 
and lie beyond the scope of this study. 
This paper is organized as follows: in section II, we will 
introduce our model based on a three-band Hamiltonian 
consisting of the Ru ^.dxy orbital and the in-plane O 2p 
orbitals. After reducing this model to an effective one- 
band model, we will analyze the effect of a rotation of the 
oxygen octahedra on it. On-site interactions are then 
treated within mean-field theory and an additional ap- 
plied magnetic field is considered. The resulting model 
is then studied in a next section. Following Metzner et 
al.,^"^ we will in section IV add a forward-scattering term 
to allow for a nematic phase and study the influence of 
the staggered spin-orbit coupling to this phase. In a last 
section we will discuss and summarize our flndings. 



II. MODEL 

A. Basic Hopping Hamiltonian 

The starting point for our model is a three-band tight- 
binding model including the in-plane 4d Ru orbital {dxy) 
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and two 2p O orbitals {px and Py) with on-site energies 
Ed and Ep = Ed — A, respectively. In this model, an elec- 
tron can hop from a d^y orbital in x (y) direction to a Pj, 
(pa;) oxygen orbital and vice versa with the hopping in- 
tegral tdpjj (tdp^)- Additionally, due to strong hybridiza- 
tion of the oxygen 2p orbitals, electrons can hop between 
neighboring oxygen orbitals. This leads to a Hamiltonian 
of the form 



n 



3b 




(1) 



pl^,p}y) are the creation operators for the 



where = (d^pt^pt 

above mentioned orbitals. Care has to be taken of the 
different signs of the hopping integrals due to the phase 
of the orbital wave-functions indicated by the tildes (see 
Fig.JI). 

To integrate out the high-energy degrees of freedom and 
thus, to reduce our model to one band, we cons truct an 
effective Hamiltonian only living at the Ru sites ,'2^'^ 



'H3b\p){p\'H3b 
Ed — Ep 



E 



n3b\p){p\n3b\p'){p'\n3b 

{Ed — Ep){Ed — Ep') 



where the sums run over all oxygen orbitals on all sites. 
This leads to a simple hopping Hamiltonian, 



(3) 



where cj^ creates an electron at Ru site i with spin s, {i, j) 
denotes nearest neighbors (nn) and next-nearest 
neighbors (nnn). The hopping integrals in this effective 
Hamiltonian then read to lowest order 



t : 



'-dp 



and t' 



f2 f 



(4) 



with tdp^ — tdpy — tdp- 

If we additionally take the rotated oxygen-octahedra 
into account, we obtain a bipartite lattice leading to a 
-\/2 X -\/2 larger unit cell separating A and B sublat- 
tices (see Fig. [s]). A consequence is that the formerly 
symmetry-forbidden hopping in x- (y-) direction via p^ 
(py) orbitals is now possible with matrix element t'^^ as 
is depicted in Fig. |4] Thus, we include spin-orbit cou- 
pling at the oxygen site, 7i'^~^^ — XL^Sz, which mixes 
the Px and py orbitals. We therefore need to change 
to eigenfunctions of the spin-orbit coupling, |±), with 
-H°~2P|±) = ±As|±), s = ±1 the spin index. 
As an example, we write the total Hamiltonian for the 
sublattice A for the x-direction. 



n 



jeA s 



i*4sP+,,-x/2s+^-^-) 



jeA s 



(5) 




FIG. 3: The Ru02 plane with rotated oxygen octahedra lead- 
ing to a doubling of the unit cell. The two distinct lattice 
sites are denoted as A and B. The inversion symmetry of 
the bonds between Ru-ions is broken leading to a staggered 
spin-dependent nearest-neighbor hopping. The next-nearest- 
neighbor hopping is still isotropic and the same for both lat- 
tice sites as is depicted by the different, equivalent hopping- 
paths (bold solid, dashed lines for A to A and dotted line for 
B to B). 



where p 



t 



creates an electron at the oxygen site 



i + x/2 in the |±) state with spin s and 



t : 



V2 



(6) 



Applying perturbation theory in the from of Eq. ^ to 
this Hamiltonian to construct an effective model, we find 
for the hopping integral from a site A to a site B in the 
positive a;-direction 



imefflB 



*\2 



A + As A - As 

y^pd ''pdJ ^2 _ ^2 

= — t + ias. 



2Xtpdt'j^d 



A2 



(7) 



Hence, we have a total Hamiltonian % — 'H^'^^ + "H*""^ 
with a nn-and nnn- hopping Hamiltonian 7^'"^ and a stag- 
gered spin-dependent hopping with the form of a stag- 
gered SOC of Rashba-type, 



ss' j^-^ o.—x,y 

+*"E E ~4.'S+-)^f.']-(8) 

jeB d=x,y 

Note that the nnn-hopping integrals, even though renor- 
malized, do not become anisotropic, which can be de- 
duced from geometrical considerations as indicated in 
Fig.|3l 

The bipartite lattice introduces a wave-vector Q = 



i^— ± as 
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FIG. 4: Possible hoppings for the three-band model in the 
case of rotated Oe octahedra. The sign change in the hopping 
integrals is due to the phase of the Wannier functions and 
holds in first order. 



(tt, tt) with which the total Hamiltonian in momentum 
space reads 



n = 



1 



EE' 



where c^^ 



-(1) 



(2) 
-k 

-*5kCrL 



-(1) 



-(2) 



My 



(9) 



' ^k+Qs^ 
COS ky — 



Here, £k = — 2t(cosfc2: + 
e^^"* + e^f' are the hopping 



- (4. 

cos ky) — At' cos kx 

energies for nn- and nnn- hopping and gk — '2a{coskx 
cos ky) is the form factor for the spin-orbit coupling. The 
prime in the k- summation indicates that the summation 
runs only over the reduced first Brillouin zone. As this 
restriction will hold for all subsequent k-sums we will 
omit the prime in the following. 

The staggered SOC hybridizes states with k and k -I- Q. 
It is now convenient to introduce Pauli matrices (t'^,t) 
in momentum space {k, k -f Q} such that we can write 
the Hamiltonian as 



^ = E 



^ks ''-kss 



ks' 



(10) 



B. Magnetization and On-site Interaction 

Before adding an on-site interaction to the Hamilto- 
nian (10), we first want to examine the effect of an ap- 



plied magnetic field. From the form of Eq. (11) we first 
see that a magnetic field in z-direction is a mere spin- 
dependent shift of the chemical potential. The response 
of the system is thus a simple polarization. However, 
for in-plane fields the staggered spin-orbit coupling in- 
troduces a coupling of homogenous magnetic fields to 
a staggered magnetization, i.e. a commensurate SDW. 
This means that the static spin-susceptibility has a com- 
ponent with wave- vector Q. 
To see this, we add a Zeeman term of the form 



n"" =g^iB[Ho-S{0) + HQ-SiQ)] 
with the spin operators 

5^(q) = jE4 



(13) 



^k+qs'^ss"^ks 



, =^5k(q), (14) 



fiB the Bohr magneton and the Lande factor g. This 
corresponds to a homogenous and a staggered magnetic 
field, in accordance with the structure of the Hamiltonian 



given in eq. (11) and can, thus, be written as 



(15) 



Here, we have introduced /ig/g ~ -^o/q/^o where Hq — 
2 • 1Q~ t / {giis)- It is now straightforward to calculate 
the magnetic response of the system to an applied field 
by using the thermodynamic relation 



d 



'O/Q/ 



-F{T,ho,hQ,N) 



(16) 



O/Q 



where (ttiq) is the homogenous magnetization pointing 
in the i-direction while (toq) corresponds to a staggered 
magnetization. Using (15) together with (10), we find 



for the case of a homogeneous field in x-direction a finite 
staggered magnetization in y-direction. 



with 



H 



(0) 



5kaL,r^ (11) 



Diagonalizing this Hamiltonian yields two (spin- 
degenerate) bands. 



+ 9i 



(12) 



where a — 1,2. The first Brillouin zone is folded back, as 
can be seen in Fig. [9^, where the Fermi surface is plotted. 
This is in accordance with the doubling of the unit cell 
introduced by the rotated oxygen octahedra. 



(™q) E EE"^(^«/5.k) 



(-l)"5k 



a=l,2B=± k 



(/»§±6W)2+52 



In the above equation, Cai.k 



,(2) 



(17) 



(-1)"^^ 



±e)y)'^+g^ are the four energy bands 



'0 4^ '') 

in a homogenous field in x-direction and 'np{^) is the 
Fermi distribution function. 



If we introduce an on-site interaction term to the 
Hamiltonian, 



(18) 



5 



which we want to treat within mean-field theory, we first 
have to choose an appropriate spin-quantization axis (in- 
dicated by the superscript a). For simpUcity, we only 
consider the two cases of an applied magnetic field in z- 
and in a:- direct ion. 



1. Field applied in z-direction 

Since a field perpendicular to the plane does not couple 
to any staggered magnetization, the first case is straight- 
forward. The quantization axis is the 2;-axis and we write 



(19) 



Since we do not expect large fluctuations in the charge 
density, the first term is a constant and, thus, the inter- 
action can be written as 



(20) 



with 5*1 = (ui-^ — ni^)/2. Applying mean-field theory to 
this expression and changing to momentum space yields 



= -2UK'P S^{Q) + UN{M 



z\2 



(21) 



with Af^ = (S'f ) independent of site i. Therefore, this 
leads to an additional term in ( 11 ) 



The total Hamiltonian now reads 



with the effective magnetic field Hq — 
plicity. This Hamiltonian has four eigenenergies 



(22) 



(23) 



U for sim- 



Co 



J2) 



+ sh^, + {^l)ygl + {ei'^y-f, (24) 



and via the grand-canonical potential per lattice site, 

-T E Elog[l + ^^P(~^"ks/r)] + C/(Af-)2, (25) 



a=l,2 ks 



the self-consistency equations can be derived, 

^ ^ Xl'^^^^^ks), (26) 

Q=l,2 ks 

2^ E E^"J=^(^"ks). (27) 



a=l,2 ks 





FIG. 5: Real space schematic of the magnetic order for the 
case of an applied field in the ly-plane: Due to the staggered 
spin-orbit coupling, the magnetization is also staggered with 
respect to the sublattice sites A and B with order parameter 
M and 6 for the total moment and canting angle, respectively. 




1.32 1.33 1.34 1.35 0.5 1 1.5 2 
n p{E) 

FIG. 6: Left Figure: Critical interaction strength for the fer- 
romagnetic (solid line) and the SDW instability (dashed line) . 
The crosses denote the choices for the interaction strength 
used for the two densities in Fig. |14[ Right Figure: Density 
of states in the absence of an external field showing the van 
Hove singularity at /i = 4t' = 1.44t (dashed line). 



2. Field applied in x-direction 

The case of a field in x-direction is slightly more in- 
volved since we expect an additional response in form 
of a staggered magnetization in y-direction leading to a 
canted magnetization in x-direction. Therefore, we de- 
fine an angle 9 for the canting angle and denote with M 
the total moment (see Figure [s]). The quantization axis, 
a and h for the spin on lattice sites A and B, respec- 



tively, in Eq. ( 18 1 should therefore be perpendicular to z 



and tilted away from the x-direction by the angle 9. We 
start from (l20|) and decouple to find 



= UNM^ - 2t/ A/ E ^ + '^i^ ^) 

-2UMY^{S^ cos 9 - 5*^ sin 9). (28) 

Changing again to momentum space, this results in an 
additional term in the total Hamiltonian 

= -(C/Af cos 6l)cr="r° - (C/Af sin 6i)cr«r^ -I- UNM"^ 

(29) 



6 



which then becomes 



Hks.' = n'^l + K<j^T^ + m (30) 



with — —UM sin 6 and 



hl^hl- UM cos 0. 



(31) 



The eigenvalues read 



6 



2±,k 



J2) 



.(2) 



(m^'i.gO' + l/iSiek^)'-^ 
(m?'±5k)2 + (/ig±e^'^)2-^ 



Again, the self-consistency equations can be deduced 
from the grand-canonical potential and are 



N 

C 

UM 



=± k 



a p=± k 



a 0=± k 



dM 



with 

561 
dM 



(-1) 



= (-1) 



,TgkCosg-(/tg/2T4'^)sing 
(m^/±gk)2-t-(/ig±ei'^)2 ' 
„ (/ig/2± e^^^ ) cos 61 - f/M±gk sin ^ 
(m?'±5k)2 + (/ig±eW)2 



(32) 
(33) 
(34) 

(35) 
.(36) 



From Eqs. (33 1 and (35) it follows trivially that there 



is only a canted magnetization if there is a finite SOC. 




2 3 4 
H^Ho 



FIG. 7: Magnetization for n — 1.336 for an applied field in 
2-direction for a temperature of T = 5 • 10~*i. Inset: the free 
energy / for H = He as a function of magnetization between 
= 0.005 and 0.01 to emphasize the first-order character 
of the metamagnetic transition. 




FIG. 8: Uniform magnetization Af^ (solid line) and staggered 
magnetization AI^ (dashed line) for n = 1.336 for an applied 
field in a;-direction for a temperature of T = 5 • 10~''t and 
a = 0.05t. For this temperature, the metamagnetic transition 
is clearly of first order. 



III. RESULTS 



Field in ^-direction 



Since a metamagnetic transition is expected to occur 
close to magnetic instabilities, we investigate their oc- 
currence from Eq. (27). The linearized self-consistent 



equation yields the condition 



U 
N 



f^4rcosh2(^,ks/2r) 



= 



(37) 



for the occurrence of a ferromagnetic instability, the fa- 
miliar Stoner criterion. This is not surprising since, apart 
from the folding of the Brillouin zone, the staggered 
spin-orbit coupling in the case of a magnetic field in z- 
direction only leads to a renormalization of the nearest- 
neighbor hopping (c.f. Eq. (24)). The critical interac- 



a function of the electron density, n, is shown in the left 
part of Fig. |6] (solid line). It drops significantly close to 
a density of UyH ~ 1.35. This corresponds to fiyH = 4i' 
where the Fermi surface hits the van Hove points located 
at (±7r, 0) and (0, ±7r) thus leading to a diverging density 
of states. This divergence is shown in the right part of 
Fig. |6j where the density of states in zero field is plotted. 
For a further analysis of the metamagnetic transition 
with applied field in z-direction, we fix the density of 
electrons slightly below UyH, n = 1.336, and choose an 
interaction strength U close to the critical one obtained 
from the linearized self-consistency equation (see Fig.|6|. 
Here and in following numerical calculations, we keep the 
spin-orbit coupling strength at a = 0.05t. The magneti- 
zation curve obtained is shown in Fig. [Tj To emphasize 
the first-order nature of the transition at T = 5 • lO^^t, 
the inset shows the free energy at the critical magnetic 
field as a function of magnetization M^ for values be- 
tween 0.005 and 0.01. 



tion strength for a ferromagnetic instability to occur as 
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Field in a;-direction 



linearizing the self-consistency equation (34), leading to 



For the case of the field applied in x-direction, we can 
again explore the occurrence of a magnetic instability by 



U 
N 



[EE 

all k 



(.gk sin 



(1) / 
ejj. cos ( 



EE'^^(^"/5'k)(' 

aP k 



-^xa (gk cos6l-£^^^sin6>) 

[^2 + (,(l))2]3/2 



4rcosh^(C„^,u/2r) gl + ie'i^'>r 



0. (38) 



To further analyze this, it is useful to change to 
different order parameters, from {M, 6) to (M^ — 
Mcos0,My = Msin6') with the uniform magne- 
tization in a;-direction and M'^ the staggered component 
in j/-direction, respectively. Writing the self-consistency 
equations in these new parameters, 



My 



^EE^^^^^^/s.k) 



a,l3 k 



(m2/±gk)' + (^g±ek')' 



1 Y-V- , i-ir{±gk-UMV) 

a,0k J('^^/±gk)2 + (/ig±eW)2 



and linearizing this system of equations. 



My 



leads to two different possible magnetic instabilities. 







1 



U 



1 



iV^4Tcosh2(ea/3,k/2r)' 



^ v/^^+gv ir'^?"^^"'''^ 

TV ^2 cos fc. 



ctk 



- COS fe, 



(39) 

h 

(40) 
(41) 



The first one is the same as Eq. (37) and corresponds 
to a ferromagnetic instability. The second equation 
corresponds to a SDW instability occurring due to the 
near nesting of the Fermi surfaces. Note, however, that 
in both cases, the magnetization will have a uniform as 
well as a staggered component. 

The solutions of these equations as functions of the 
electron filling, n, are plotted in Fig. [6j We see that the 
critical interaction strength for the spin-density wave 
instability is generally below the one for a ferromagnetic 
instability, but shoots up when approaching the critical 
filling. 



We can now interpret the linearized equation ( 38 ) as 



having a ferromagnetic and a spin-density wave contri- 
bution. Proximity to a SDW instability additionally 
lowers the critical interaction strength. Therefore, the 
metamagnetic transition could occur at a lower field 
then in the z-direction case, especially close to a SDW 
instability. 



The magnetization due to a magnetic field applied in 
plane, as well as the amplitude of the staggered magneti- 
zation, are plotted in Fig. [Hj Again, we find a first order 
transition for lower temperatures while the transition 
changes to a crossover upon increasing temperature. 
Note that the sign of Af^ depends on the sign of the 
spin-orbit coupling constant a. There is no degeneracy in 
the state obtained which could lead to domain formation. 



C. Comparison 

Comparing the two cases of fields applied in z- and x- 
direction in Figs. [7] and |8) respectively, we first see that 
the critical field for the latter is shifted to lower fields, 
even though the zero-field susceptibilities differ by less 
then a percent, Xo/Xo ~ 0.998. This is due to the prox- 
imity of the system to a SDW instability as was pointed 
out in the previous subsection. 

This behavior is in qualitative agreement with the exper- 
imental phase diagram (see the schematic phase diagram 
in Fig. [T]). However, the difference of the in-plane and 
out-of-plane critical field is smaller in our model calcula- 
tion than in the experiment. Our model only includes the 
staggered spin-orbit coupling entering through the oxy- 
gen displacement. Naturally, other spin-orbit coupling 
contributions, particularly from the Ru-ions, would add 
to the anisotropy through an anisotropic g-tensor, likely 
with a larger polarizability in the basal plane than along 
the z-ax\s}^ This is, however, beyond the scope of this 
study as a detailed analysis would require to include other 
bands. 

Second, a numerical study of the Gibb's free energy shows 
that the temperature T* up to which the first-order tran- 
sition persists, is higher in the case of the in-plane field. 
For our choice of parameters we find w 9 • 10^'*^ while 
T* K, 11 • 1Q~H. This anisotropy in the critical temper- 
ature is consistent with the trend in the experimental 
situation. However, it does not reproduce the quantum 
critical endpoint. Note, that the difference between T* 
and T* could not be explained simply by an anisotropic 
g-tensor. In principle, it may be possible to tune the 
model in such a way as to press the critical temperature 




FIG. 9: (a) Fermi surface for an electron density of n = 1.336 
without an applied field. Due to the rotation of the O- 
octahedra, the Brillouin zone is folded back, denoted by the 
dashed lines, (b) and (c): Fermi surface just below and above 
the critical fields in z- and s-direction, respectively. For clar- 
ity, only a small section of the BZ is shown indicated by the 
little square in (a). For the case of a field applied in plane 
(c), small gaps open close to the van Hove points. 



for out-of-plane fields, T* to zero while still having a first- 
order transition at finite temperatures for in-plane fields. 
Also fluctuation effects are likely important in this con- 
text. These features are, however, not essential to our 
discussion. 

An additional important difference between in-plane and 
out-of-plane fields can be seen in Fig. [9] where the Fermi 
surfaces for both cases for fields below and above is 
shown: we see that the system undergoes a metamag- 
netic transition to prevent the majority-spin band from 
touching the van Hove points. In (c), we additionally see 
that the induced spin-density wave opens small gaps at 
the Fermi level close to the van Hove points. This has 
important consequences for the appearance of a nematic 
phase as we will see in the next section. 



FIG. 10: Mean-field results for the magnetization and the 
nematic order parameter r\ for a field applied in ^-direction 
for a temperature of T = 5 • 10~*t showing an intermediate 
nematic phase bounded by two first-order transitions. Here, 
g = 0.33 and n = 1.336. 

IV. NEMATIC INSTABILITY 

In this section, we explore the occurrence of a nematic 
phase in our model for the two cases of a magnetic field 
applied in z- and x-direction, respectively. For t his pu r- 
pose, we introduce an additional interaction ternpSEU 

^" = ^EE/kk'nk.nk'.' (42) 

kk' ss' 

with a coupling function /kk' only contributing for zero 
momentum transfer, i.e. for the forward scattering, 
which is the relevant interaction for a nematic phase to 
occmW^ We then separate the coupling function 

/kk' = ffrfkc^k' (43) 

and choose a A.j.i_yi symmetric form for the form factors, 
dk = cos /ca;— cos ky. This term can then lead to a nematic 
phase, reducing the symmetry from C4 to C2. 
Introducing again a mean-field decoupling which is spin- 
independent we write for this interaction 

H" = ;^77dknks-^ry' (44) 

ks ^ 

with 

E'^k'(nk'). (45) 

k' 

Since = — rfk+Qj but is isotropic in spin-space, we can 
deal with it by replacing 

(46) 

while all the above formulae still hold with the additional 
self-consistency equation 
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FIG. 11: Phase diagram for a magnetic field applied in z- 
direction. While for low temperatures, the two consecutive 
transitions are of first-order (solid line), they become second- 
order before the nematic phase disappears completely (dashed 
line). Above this temperature, a metamagnetic crossover can 
still be seen. 



for the z-direction case and 



'7 = -^I]I]"F(W,k) 

a,l3 k 



(-irdk(ef^±/i§) 



/{rhy±gur + ih-o±el'>)^ 

(48) 

for the x-direction case, respectively. 
For sufficiently strong we find a magnetization curve 
for fields applied in 2;-direction as is shown in Fig.[lO) The 
two jumps in the magnetization border an intermediate 
phase with a finite value of the nematic order parame- 
ter rj. The instability in this case is again driven mainly 
by electrons whose momenta lie close to the van Hove 
points. To obtain an intermediate phase before a single 
metamagnetic jump removes all such electrons from the 
Fermi surface, a critical scattering strength gc is neces- 
sary. Above that, a nematic phase is entered at some 
magnetic field Hd and left again at Hc2- 



The T-iJ-diagram shown in Fig 11 shows first-order tran- 
sitions for low temperatures up to T « O.OOli, second- 
order transitions for higher temperatures until at T w 
0.0016t the nematic phase disappears completely to make 
way for a metamagnetic cross-over (not shown). This 
beha vior h as already been observed in similar calcula- 
tions.l™ 

For the case of a field applied in a:;-direction, a very sim- 
ilar behavior is observed, however, with one important 
difference: Since the induced spin-density wave already 
removes some weight from the Fermi surface close to the 
van Hove points, a larger forward scattering strength 
is required for the occurrence of a nematic phase, i.e. 
9c > 9c- Therefore, there is a range of g, where there 
exists already a nematic phase for fields in z-direction, 
but only one metamagnetic jump is observed for in-plane 
fields. This result is summarized in the phase diagram 
in Fig. [T2j As a function of the SOC strength a and 
the forward scattering strength 5, we find three regions. 
In addition to the two obvious ones, where there is ei- 




0.3 0.325 

9 [t] 



0.35 



FIG. 12: Critical forward scattering strength for the cases of 
an applied field in z- and x-direction, respectively. We can 
distinguish three different regions. I: The forward scattering 
strength is too week to enter a nematic phase, no matter in 
what direction the magnetic field is applied. Only a single 
metamagnetic transition occurs. II: While there is a nematic 
phase for fields applied in z-direction, no such phase occurs 
for fields in plane. This region corresponds to the case found 
in Sr3Ru207. Ill: the forward scattering strength is strong 
enough such that a nematic phase will occur for fields in any 
direction. Inset: schematic of the phase diagram showing that 
the nematic phase only occurs for fields close to the z-axis. 



ther no intermediate phase (region I) at all or one for 
fields applied in any direction (III), there is now a new 
region with a nematic phase only for fields applied in z- 
direction (II). Obviously, this region corresponds to the 
case of Sr3Ru207. 



V. DISCUSSION AND CONCLUSIONS 

The range of the forward scattering strength for which 
in our calculation a nematic phase is only observed for 
fields in z-direction is not very large. Note , however, that 
the nematic phase only appears in a very narrow region 
((ffc2 - Hci)/Hai < and thus, g is only slightly 

bigger than g^ 



This is illustrated in Fig. 13 where the 
dependence of the critical fields on the forward scatter- 
ing strength is shown. The width of the nematic phase 
grows rather rapidly with increasing g. Therefore, the 
actual size of the forward scattering strength might well 
lie in region II of Fig. [T2j 

Obviously, not only the strength of the SOC, but also 
the nesting properties of the spin-polarized Fermi sur- 
faces play an important role for the appearance of the 
anisotropy effect. Nesting properties are a factor of tun- 
ing our model to the vicinity of a SDW instability. To 
analyze the impact of enhanced SDW correlations, we 
consider two different electron densities and correspond- 
ing on-site interaction strengths U with U™ — U =const 
(see crosses in the left part of Fig. 6). This allows us to 
examine the cases of two different proximities to a SDW 
instability with comparable strengths of FM correlations. 
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FIG. 13: Phase diagram for a magnetic field applied in z- 
direction depending on the forward scattering strength g. Be- 
low a critical strength gl , no nematic phase is entered, but the 
system undergoes a single metamagnetic transition connected 
to the proximity to the ferromagnetic instability as discussed 
above. Above g^, the system enters a nematic phase whose 
region grows with increasing g. Here, T = 5 ■ W~'^t. 



One important finding is that the stronger the SDW cor- 
relations the more pronounced the anisotropy effect and 
thus, the smaller the ratio H^/H^ becomes. This is de- 
picted in Fig. 14. In Fig. 6 we also show that the ferro- 
magnetic and SDW instabilities can be competing for the 
chosen parameter range. We fix our model parameters in 
a way to avoid the occurrence of a staggered magnetic 
moment for any value of the magnetic field along the 
z-axis, while the staggered moment is field induced for 
in-plane fields. 

We should also comment on the strength of the SOC a 
that we expect for this system. To get an estimate of the 
on-site SOC strength A for p-electrons on the oxygen, we 
take the vacuum value, A ~ lOmeV. A very crude 
estimate of the staggered SOC coupling from Eq. ([7| then 
yields 



a 
7 



2Ai;^ 
A U, 



2A 



(49) 



Taking estimates for A k, 1.5eV we find that for a a value 
on order of a percent of t seems reasonable .I^Si Comparing 



Figs. 12 and [T3j this would allow for a nematic phase 
with a width of (i?c2 — Hci)/Hci « 2%, in agreement 
with experiment. For a more reliable estimate of a, DFT 
calculations should be performed. 

Finally, some remarks to the nematic phase are in order. 
As was already mentioned, the nematic phase introduced 
here is the same as discussed by other authors.l^^lti^ As 
shown by these authors, the nematic phase could account 
for several experimentally observed phenomena, like the 
anomalous resistivity, or the non-Fermi-liquid behavior 
of the susceptibility and the specific heat coefficient. 

To conclude, we showed that the rotated oxygen oc- 
tahedra lead to a staggered hopping that can be de- 



scribed with the help of a (staggered) spin-orbit coupling 
of Rashba-type. This introduces an anisotropy of the 
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FIG. 14: Anisotropy of the critical field in a;-direction com- 
pared to the value for the z-direction for the densities n = 
1.336 (solid line) and n = 1.339 (dashed line), respectively. 
The value for the on-site interaction is chosen as indicated in 
Fig. |6] The dependence on the strength of the spin-orbit cou- 
pling is more pronounced the closer the system is to a SDW 
instability. 



response to a magnetic field, namely an induced spin- 
density wave for the case of in-plane fields. This stag- 
gered magnetization could be observed in neutron scat- 
tering experiments. To our knowledge this kind of ex- 
periment has not been performed so far. The additional 
magnetization has, first, the effect that the critical field 
for a metamagnetic transition is shifted to lower values 
for in-plane fields. Also, the critical temperature T* up 
to which the transition is first order is higher for fields 
in the xy-plane. Last and most important, the spin- 
density wave opens gaps at the Fermi level that lead to an 
anisotropy for the appearance of a nematic phase. Addi- 
tionally considering spin-orbit coupling effects of the Ru 
orbital would account for the full anisotropy of He {g- 
tensor anisotropy). Therefore, the present work allows 
for a picture that is qualitatively consistent with experi- 
mental observations including the anisotropics in H^, T* 
and the appearance of a nematic phase. 
Note added: While preparing for submission we noticed 
the recent papei^^ which studies the inffuence of spin- 
orbit coupling and the doubling of the unit cell on the 
nematic phase of Sr3Ru207. 
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